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Abstract: Observational evidence indicating that the expansion of the universe is accel- 
erating has surprised cosmologists in recent years. Cosmological models have sought to 
explain this acceleration by incorporating 'dark energy', of which the traditional cosmo- 
logical constant is just one possible candidate. Several cosmological models involving an 
evolving equation of state of the dark energy have been proposed, as well as possible 
energy exchange to other components, such as dark matter. This paper summarises the 
forms of the most prominent models and discusses their implications for cosmology and as- 
trophysics. Finally, this paper examines the current and future observational constraints 
on the nature of dark energy. 
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1 Introduction 

The observed acceleration of cosmic expansion is a 
landmark discovery in modern Cosmology. We don't 
know why the universe is accelerating, but it is reason- 
ably clear that the explanation will require new fun- 
damental physics — beyond our current understanding 
of particle physics, gravitation and the quantum vac- 
uum. Cosmologists are now in the fortunate position 
of being able to ask the questions that will provide in- 
sight into ultimate physics beyond the reach of particle 
accelerators, as well as understanding the history and 
fate of the universe. 

The present best fit cosmological model, known 
as the concordance model, combines data from many 
complementary sources incl uding the WM AP and other 
observations of the CMB llSpergel et al. ' 2003), l arge 
scale structure surveys such as the 2DFGRS ( H a wkins et 
I2003D and the SPSS llAbazaiian et alJl2003l) and su - 
pernovae data jKnop et alJl2003l: iRiess et alJl2004D . 
The accelerated expansion of the universe is modelled 
via the 'cosmological constant' (A), an entity that has 
been introduced, removed and re-introduced several 
times since Einstein's original introduction into his 
field equations. The cosmological constant has some 
special properties that make it a natural choice for in- 
clusion into our models, however there are many plau- 
sible alternatives also permitted by current data. The 
purpose of this paper is to review cosmological model 
basics and then explore some general unheralded im- 
plications of the many cosmological models that utilise 
a generic 'dark energy' rather than assuming a priori 
a cosmological constant. 

While there is good agreement between current ob- 
servational data and cosmological models incorporat- 
ing a non-zero cosmological constant, its introduction 
has not been without problems. The first is known as 
the coincidence problem. The energy density of matter 
decreases as the universe expands, proportional to the 



cube of the scale size of the universe. However, the en- 
ergy density associated with the cosmological constant 
remains constant as the universe expands. Thus, that 
we should be observing at an epoch when the energy 
density of matter and the cosmological constant are of 
the same order of magnitude seems very unlikely. 

The second problem is the infamous cosmologi- 
cal constant problem, discussed in detail in IWeinberd 
(198^) and |£arroll (2001). In short, the fact that 
the energy density of the cosmological constant is un- 
changed by the expansion of the universe suggests that 
it can be identified as a property of spacetime itself - a 
vacuum energy density. Such a zero-point energy can 
be calculated from quantum field theory. However the 
theory predicts that if the vacuum energy is not zero, 
then it should be a value that is 120 orders of magni- 
tude greater than the energy density required for the 
al] cosmological constant. 

The problems with the cosmological constant, as 
well as the relative lack of observational constraints, 
have lead to a flurry of alternative explanations in re- 
cent years. A component that causes the expansion of 
the universe to accelerate is referred to as dark energy, 
with a cosmological constant being just one possibility. 
Another popular candidate is a primordial scalar field, 
or 'quintessence'. This is the generic name for a time 
varying, spatially inhomogeneous component with neg- 
ative pressure. The prime example is a scalar field Q 
slowly rolling down a self- interaction potential V(Q) 
(Ratrafc Peebles 1988. : Steinhardt fc Caldwell 199j ; 
ISteinhardd l200a) . Unlike the cosmological constant, 
which is ascribed to vacuum energy, the quintessence 
field has no expected value and the 120 orders of mag- 
nitude problem can disappear. However, it is a much 
more ad hoc approach as there is no 'natural' reason 
(besides possible connections with early universe infla- 
tion) to postulate the existence of a quintessence field, 
unlike the more comfortable vacuum energy interpre- 
tation of the cosmological constant. Many and varied 
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Figure 1: Cosmological evolution for a constant 
equation of state: the top two panels show the en- 
ergy parameter for the dark energy component Oa, 
the middle two panels are for matter fim and the 
bottom two plot the scale factor. The left panels 
show evolution with respect to time (in units of 
1/Hq), while the right panels are with respect to 
redshift. The bottom right panel plots an identity 
(equation (jSJ) and is shown for consistency. 
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Figure 2: Distance measurements for a constant 
equation of state: the top two panels show the ef- 
fective magnitude as defined by equation (|18|l . us- 
ing the angular diameter distance in units of Mpc 
using an arbitrary v alue of h=0 72. Th e gold and 
silver data sets as in iRiess et all (1200411 are shown 
in corresponding colours. The bottom two show 
the angular diameter distance in units of 1/Hq. 
The left side panels show evolution with respect 
to time (in units of I/Hq), while the right panels 
are with respect to redshift. 



quintessence models have been proposed l^ ggg^rong^ 
l2005l : IChimento et alj2003l:lBarenboim fc Lvkkenl2005fl . 

The evolution of the properties of the scalar field may 
solve the coincidence problem, as described in section 
15.21 There are also alternative models of gravity that 
seek to explain the obser ved data [for instance see 
llCarroll et alJ |2004 l2005t) and references therein] as 
well as explanations other than accelerated expansion. 
If distance supernovae were dimmed by some kind of 
'grey dust' then this would give the false impression 
cosmic acceleration, however see ()R,iess et alJ20oJl for 
a discussion of the problems with the grey dust model. 
More exotic dimming pro cesses such as ph oton /axion 
mixing in magnetic fields ijCsaki et all2003l have been 
proposed, though calcualtions show that this cannot 
alone account for the observed acceleration. Observa- 
tional data are currently insufficient to cull the field of 
alternatives. The solution of this critical problem in 
modern cosmology will require a great deal of obser- 
vational efi'ort. 

This paper examines the general form of several 
of the currently proposed models, including those in- 
volving energy transfer between dark energy compo- 
nents, and examines their infiuence on the past and 
future expansion of the Universe. Section |21 gives a 
general background on cosmology, introducing the no- 
tation and important equations as well as outlining the 
current observational evidence. Section |3 introduces 
the equation of state and examines models where it 
is allowed to evolve with time. Section |1] introduces 
models with an interaction between dark energy den- 



sity and matter density. Section |3 discusses the im- 
plications of the models presented in sections |3| and 

H 



2 Background 

Modern theoretical cosmology is built on two pillars: 
the cosmological principle and general relativity. The 
cosmological principle states that we do not occupy a 
special position in the universe. This allows our obser- 
vation that the universe is isotropic to be extrapolated 
to the global property of homogeneity. Homogeneity 
dramatically simplifies our cosmological theories, as it 
means that any parameter that describes the universe 
as a whole can only depend on time. 

From the assumptions of homogeneity and isotropy, 
a metric can be derived that tells us how to mea- 
sure distance and time in the universe. Thi s metric is 
known as the Robertson-Walker metric (see IWeinberd 
(|l972fl for more details), and has the line element: 

ds^ = c'&t^ — 

B? {t) {dx^ + S'i [AO'^ + sm ed<f')) (1) 
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Figure 3: Cosmological evolution for an equation 
of state linearly varying in redshift: the layout is 
the same as figure ^ Note that this parameterisa- 
tion breaks down for z^ \. 




Figure 4: Cosmological evolution for an equation 
of state linearly varying in redshift: the layout is 
the same as figure |5] Note that this parameterisa- 
tion breaks down for z > 1. 



where c is the speed of light (hereafter c = 1) and 

{sinx closed universe (fc — +1) 
X flat universe (fc = 0) (2) 

sinhx open universe (k = -1) 

where (x, ^, </>) are the spherical comoving coordinates, 
which for an object moving with the Hubble flow do 
not change as the universe expands. 7? is the scale 
factor of the universe, which is the key prediction of 
any cosmological model. It contains information about 
evolution and fate of the universe, and is simply related 
to the observable redshift (2) by: 



Ro 



1 + z 



(3) 



The subscript refers to the present epoch. We also 
define the Hubble parameter H, which measures the 
rate of expansion of the universe: 



H 



R 
R 



Differentiating ||3J gives: 



H 



dz 



1 + z dt 



(4) 



(5) 



In the curved, expanding spacetime of the Robertson- 
Walker metric, no t all methods of measuring di stance 
are the same (see iLinderl II1997D for details and lHoed 
(^99) for a summary). The comoving radial distance 
(rp) between the origin and (x, 0, 0) is defined as a si- 
multaneous {dt = 0) radial measurement {d9 = d(j) = 
0) at time to and is given by: 



(6) 



Light travels along null geodesies, defined by ds — 0. 
Since we have placed ourselves at the centre of the 



coordinate system, light moves radially {dO = d(/) = 0). 
Noting that all the quantities are positive, the RW 
metric reduces to: 



cdt = -R(t)dx 



X(iem) = C 



dt' 



(7) 



Note that this is not simply multiplying the speed of 
light by the light travel time (to — tem). To convert 
equation Q in terms of the observable redshift, we 
use 

where we have used the fact that to corresponds to 
2 = and tern to z. 

Proper transverse distance (rm) and angular diam- 
eter distance (r^) both relate the transverse proper size 
I of an object at x to the angular size {A9). For Vm, 
all measurements are done simultaneously (dt = 0) at 
time to giving: 



= RoSk{x) 



(9) 



For Td, measurements are done using photons which 
travel along geodesies (ds — 0). Thus, angular diam- 
eter distance to an object which emits light at time 
(tern) as measured by an observer at to is the proper 
transverse distance at tern'- 



Td = R{tem)Skix) 



(10) 



Luminosity distance (ri) to a source of intrinsic 
luminosity L from an observer who measures its fiux 
S is defined by: L = 47rr-f S. The distances rm, and 
ri are related by: 



n = (1 + z)rm ^ {1+ Z) Td 



(11) 
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The field equations of general relativity allow us 
to relate the metric of the universe Q to its energy 
content. The result is the Friedmann equations: 



P 

R 
R 



BttG k 

-3H(p + p) 
4nG 



(p + 3p) 



(12a) 
(12b) 

(12c) 



where an overdot refers to differentiation with respect 
to time; G is Newton's gravitational constant; p is the 
total energy density; p is the pressure; H is the Hubble 
parameter; and c = 1 for convenience. 

Equation H12a^ is known as the expansion equa- 
tion, Ill2bll is the adiabatic equation and equation 112cl 
is the acceleration equation. Any of these equations 
can be derived using the other two, although all three 
equations are produced separately from the Einstein 
field equations of general relativity (and are related 
through the Bianchi identity). 

As is consistent w ith CMB data, including WMAP 
llSoergel et al.ll2003[) . we assume throughout that the 
universe is flat (fc = 0). We consider the universe 
to contain a number of components (labelled i), each 
with a corresponding pressure pi and density pi that 
contribute to the total p = X^iPii P ~ X^i Pi- Defining 
the critical density to be: 



Pcrit(2;) 



SttG 



(13) 



the energy components can then be described relative 

to pcrit: 



(14) 



pcrit (z) 

where fli is the dimensionless density parameter for 
component i. Then, by putting fc = into l|12a|l the 
sum over all the density parameters is unity: 



(15) 



Details of how the density parameters ili{z) were 
calculated in this paper are contained in appendix IXI 

As well as energy density, pressure also appears 
in the Friedmann equations. We commonly relate the 
pressure of a component to its energy density using an 
Equation of State (EOS): 



EL 

pi 



(16) 



If the components of the universe are non-interacting 
(they do not exchange energy), then they will each 
satisfy their own adiabatic equation (I12bt : 



Pi = -3Hpi{l 



(17) 



For most of the familiar components of the universe, Wi 
is constant: ordinary, non-relativistic matter is essen- 
tially pressureless and has lOm = 0, radiation h as Wi = 
1/3, a cosmological constant has iua ~ — 1 l|Linderl 
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Figure 5: Cosmological evolution for the linear pa- 
rameterisation of the equation of state detailed in 
section l5.2.2l the layout is the same as figure ^ 



A cosmological model, capable of specifying R{t) 
for all times, is described by specifying: the compo- 
nents of the universe, their equations of state Wi, their 
densities today Qi (0) and the present value of the Hub- 
ble parameter, Hq. Experimental uncertainty in Hq is 
often written in terms of the dimensionless h, defined 
by _ffo — 100/iMpc~^kms~^. The dependence of quan- 
tities on Ho can then be made explicit, e.g. Qih^ . The 
value of h was measured by the Hubble Key P roject 
team to he h = 0.72 ± 8 llFreedman et al.ll200lll . 



2.1 Observational Evidence 

Several independent lines of evidence lead to the con- 
clusion that the expansion of the universe is accelerat- 
ing. Recent observations of the CMB from the WMAP 
satellite have confirmed that the universe appears to be 
flat, with n = 1.02 ± 0.02 fSocrecl ct al. 2003). The 
power spectrum of galax y distributions on th e other 
hand points to ^ 0.3 ijPercival et al.ll2002^ . A key 
third measurement is the magnitude-redshift relation 
of supernovae type la (hereafter SNla), see Riess et ^ 
(2004) for a compilation of recent data. These mea- 
surements give strong evidence that the expansion of 
the universe is accelerating, and has been accelerating 
since z ~ 1. The concordance model accounts for the 
missing energy density needed to make the universe 
flat via a cosmological constant {w = —1), however 
this is by no means the only permitted form of dark 
energy. A new generation of cosmological experiments, 
including the proposed SuperN ova Acelleration Probe 
(SNAP) jAldering et al.ll2003) are expected to mea- 
sure more supernovae at higher redshifts. This new 
data will be crucial in discriminating between cosmo- 
logical models. 

SNla are very bright standard candles, making 
them excellent cosmological probes. The distance mod- 
ulus, once the appropriate K-correction has been ap- 
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Figure 6: Cosmological evolution for the linear pa- 
rameterisation of the equation of state of detailed 
in section 15. 2. 21 the layout is the same as figure |5| 
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Figure 7: Cosmological evolution for a constant 
interaction between matter and dark energy: the 
layout is the same as figure ^ 



plied is given by: 



M = 25 + 51og(i/on) 



(18) 



where ri is defined in section |5| and is presented here 
in units of Mpc. AU three distances previously de- 
fined are model dependent and hence the measure- 
ment of the magnitude and redshift of SNla can be 
used to discriminate between models. The methods 
used to compute distances are described in appendix 
1X1 The results are presented in the even numbered 
figures i n this paper . In th ese figures SNla data com- 
piled in lRiess et al.l l|2004fl have been added in order 
to illustrate how similar all the models are in the re- 

[ fon z < 0.5 where the data sit. The SNAP probe 
Aldering et al.ll2004fl is anticipated to reach z ~ 1.5, 
where the models can be seen to diverge much more 
than for the current data set (see even numbered fig- 
ures), illustrating the much greater constraints on the 
models that will be possible. 

The data described above lead to a cosmological 
model where 

fc = (The universe is flat) 



h ~ 0.7 iln 



: 0.3 ^Ao ^ 0.7 



Note however that this model makes the assumption 
that the dark energy is a true cosmological constant. 
This assumption is permitted by existing data, but not 
strongly favored over other possibilities. Despite this 
untested assumption this model is used widely. In this 
paper the universe is assumed to be fiat with QmO = 0.3 
and f^xo ~ 0.7 as in the concordance model. Note that 
the subscript X is used to denote dark energy, leaving A 
to denote the specific case of a cosmological constant. 
Ho is set to unity, which is equivalent to measuring 
time in units of H^^. 



3 Equation of State of Dark 
Energy 

In this section the effect of different forms of the dark 
energy EOS (as defined by CSJ) on the expansion his- 
tory and distance measurements are examined. We 
take a phenomenological view of the EOS, rather than 
speculating on the physical processes behind any par- 
ticular form. This is a common approach given the 
mysterious nature of dark energy and current data be- 
in g unable to constra in com plex models. It was shown 
in lLinde~ fc JcnkinJ i2003t) that any term other than 
matter in H12a|l can be modelled by an effective equa- 
tion of state of dark energy, regardless of the physical 
origin of the term. If we write (I12al as 



(19) 



with Ho set to unity, then by comparison to the case 
of dark energy with some equation of state w 



H^ = finii?" -h (1 - Q.u.)R 



-3(1-1- lu) 



and 



^ _ ld(\n&H'^) 



(20) 



(21) 



3 d(ln_R) 

For any arbitrary mechanism that causes there 
is a corresponding function w (not necessarily con- 
stant, as in H2()^ 'l that can be calculated. This allows a 
wide variety of different models and mechanis ms to be 
compa red within a single parameter space fsee iLinded 
II2004D '). 

3.1 Constant Equation of State 

Current data are ins ufficient to constrain w x to more 
than one parameter llLinder fc MiaueJl2004D : in other 
word s, a constant wx. Using C MB, LSS and SNla 
data, IWang fc Mukheried i2004ll found the constraint 
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Figure 8: Cosmological evolution for a constant 
interaction between matter and dark energy: the 
layout is the same as figure HI 
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Figure 9: Cosmological evolution for a decaying 
dark energy model: the layout is the same as figure 
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on a constant dark energy EOS to be —1.24 < uix < 
—0.74. We present a wider range of possibilities tlian 
this, in order to exaggerate the effects for clarity. When 
Wi is constant, equation 11711 can be easily integrated 
to give Pi = /5i(z) using equation Q: 



(22) 



Equation 12211 gives the familiar results that, as we look 
back into the universe, pm cx (1 + zf" , pr oc (1 + z)*, 
whilst PA remains constant as the universe expands. 
From we can see that as the universe ages, the 
component with the most negative w will come to dom- 
inate the total energy of the universe. Consider equa- 
tion H12c|l in the case where the component with the 
most negative w has come to dominate the total energy 
of the universe. If ui < — |, then R will be positive and 
the universe will accelerate. We can therefore divide 
components into accelerating (w < — |) and deceler- 
ating (it;>i). Ifiu = — i then the expansion coasts: 
a =constant. 

Thus, if the universe has any accelerating compo- 
nents, the universe will eventually begin to accelerate, 
independent of its geometry ^ . Only if all the compo- 
nents are decelerating will the universe decelerate, and 
its fate be decided by its geometry. 

Things get interesting when w < — 1. If we con- 
sider an equivalent form of equation 1221 : 



Pi(i?) = piO 



-3{l+ii.i) 



we see that w < —\ implies that the energy den- 
sity of the component will increase as the universe 



^The exception is where the energy density in the de- 
celerating components is high enough that the universe be- 
gins to contract before the accelerating components come 
to dominate the total energy density. 



expands. Such a component has been dubbed "phan- 
tom energy" and its consequences have been studied in 
ICaldwell. Kamionkowski. fc Weinberel (;2QQ3.1 . In short, 
phantom energy causes the scale size of the universe 
to approach infinity in finite time. The resulting "i? = 
CX3" singularity has been dubbed "the big rip". Before 
the end, the expansion of the universe will overpower 
all other forces - first gravity, as the Milky Way, the 
Solar System and finally the earth are pulled apart, fol- 
lowed by weak, EM and strong nuclear forces as matter 
is torn into its components. 

In figures and |5| we consider models identical to 
the concordance model (matter -I- dark energy) except 
that we allow the equation of state for dark energy to 
vary. Figure shows the expansion history and den- 
sity parameters while |5]shows the measurable distance 
properties with the current SNla data overlaid. 

The pink model [w — —4/3) is an example of a 
phantom energy model. The "big rip" behaviour as 
time increases can be seen in the bottom left panel of 
figure^ The light blue model {w — —1) is the concor- 
dance model. For the light green model {w = —0.33), 
we see from equation 112cll that as dark energy be- 
gins to dominate, the acceleration goes to zero. Thus, 
R(t) approaches a linear function as time increases. 
The red model {w — 0) is the Einstein-de-Sitter model 
(fiat, matter only), which decelerates as time increases 
but never turns around and begins to contract. The 
black model {w = 0.33) contains a radiation compo- 
nent whose energy density falls quickly as the universe 
expands. 



3.2 Evolving Equation of State 

While current data are insuffici ent to constrain wx(z) 
to anymore than one parameter llLinder fc Miauel2004D . 
next generati on data, such as the P lanck Surveyor 
^TaubeI^l2004^ and the SNAP satellite ijAldering et alJ 
1200411 will be able to constrain possible evolution of 
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Figure 10: Cosmological evolution for a decaying 
dark energy model: the layout is the same as figure 

m 



wx{z)- This section presents two different first or- 
der parameterisations of wxiz) and demonstrates the 
effects of several values of the parameters. Since an 
evolving EOS will only be considered for the dark en- 
ergy component (i = X), we omit the subscript on w 
in this section. Again using @ and Ijl2bfl 

px(^) = px(0) e'^» ^TTP^''-' (23) 

While a two parameter evolving equation of state 
w{z) could take many forms, two are most common in 
the literature. 

3.2.1 Old Linear Parameter isat ion 

The simplest first order expansion of the equation of 
state gives w{z) — wo+wiz. While obsolete (see \'A.'2.'I\i . 
this form has been used widely in the past. In this case, 
equation 123|l becomes: 

Px{z) = px(0) (1 + ^)3{l+"0— l)g3^1- (24) 

In figures 13 and 0] we show the effects of a modified 
concordance model by using a linear parameterisation 
of the dark energy equation of state. 

The first thing to note in Figurej^is the bizarre be- 
haviour of the energy densities at early time (i.e. large 
redshift). This shows the inadequacy of this simple 
parameterisation: it is inappropriate for z > 1 as w{z) 
is unbounded as 2 ^ oo (as is pxiz) for wi > 0). 

As we model into the future, further problems arise. 
For a universe that is always expanding, forward mod- 
elling is done by allowing z —1. In this case, w{z) — > 
Wo — wi. If Wo — u^i < — 1, then eventually w < —1 
and the universe will end in a "big rip" . 

3.2.2 Linear Parameterisation 

In iLindeJ (|2003t l . a parameterisation is used of the 
form: w{z) = wq + wi{l — R/ Ro) = wo + wiz/{l + z). 



Then the equation for the density becomes: 

px{z) = px{0) (1 + ^)3(i+»o+»i)e-3»ixt^ (25) 

The advantages of this parameterisation are given 
in hinder (2003), in short the old parameterisation 
clearly becomes problematic at high z, while this pa- 
rameterisation has a bounded behaviour as 2 ^ 00. 
The results of several values of the parameters in this 
model are shown in figures |3 and El 

The energy densities for the parameters shown in 
figure 13 are remarkably similar throughout the life- 
time of the universe. The past history of R{t), too, 
is very similar. However, the parameterisation is un- 
suitable for modeling the future, z —1, since w{z) 
becomes unbounded. Note by contrast that a param- 
eterisation linear in redshift is unbounded in the past 
- wher e all the data is! Models such as the "e-fold" 
case of hinder & Hutoror (2005) and the "kink" case 
of ICorasaniti et ah. (..2003.) can smoothly handle both 
the past and future. 

4 Interacting Components 

In the previous section we assumed that dark energy 
does not interact with matter or radiation. However, 
as we have no knowledge of the microphysics of dark 
energy, it is worth considering what the effects of inter- 
acting dark energ y might be. St udies of such models 
are ongoing, see ()Zimdah]l 120051 ) for a recent exam- 
ple. Much of this work has been presented in Physics 
rather than astronomical journals, as the interest cen- 
tre of the nature of the fields involved in the interac- 
tion. Here we focus the general form of such models 
and highlight the consuqunces that may be of interest 
to astronomers and cosmologists. 

The Friedmann equations do not restrict energy 
components to be non-interacting and can be easily 
modified to model an energy exchange. Equation l|17^ 
is true in the case where the component i is not ex- 
changing energy with any other components. We now 
return to equation Ill2bll and relax this assumption. 
Consider the case where there are only two compo- 
nents in the universe: pressureless matter and dark 
energy. Then ([T2bl can be split into its components 
(ISzvdlowskill200!Tft : 

{pm + 3Hpm) + {px + 3Hpx{l + wx)) ^0 (26) 

At this point, a phenomenological parameter 7 can 
be introduced to represent the interaction between dark 
energy and matter. Then: 

Pm + 3Hp^ = 7 (27) 
px + 3/fpx(l + mx) = -7 (28) 

7 has units of energy per unit time per unit volume, 
and represents the rate at which energy is transferred 
from the dark energy component to matter (7 > 0) or 
vice versa (7 < 0). 

We cannot simply integrate l|27|l and H28|l in the 
same way we did to obtain equation 12311 . To find pm. 
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px, and R{t) we need to solve l|27|l and l|28|l simulta- 
neously with the expansion equation I12at . This rep- 
resents a coupled set of (potentially) non-linear ODE's 
which are solved numerically (see appendix EJ. 

Note that in this section we take wx = — 1 for all 
the models, though this is not a cosmological constant. 
We do this to isolate the effects of the interaction on 
the evolution. Energy transfer between matter and 
other forms of dark energy could also be computed, 
however we are interested in examining the broad con- 
sequences of energy transfer, which will be qualita- 
tively similar for other forms of dark energy. Future 
work will generalize this; see also iLindeil lj200firi . 

4.1 Constant Interaction 

The first form of interaction considered was a constant 
transfer of energy, 7 = constant. We considered only 
the case where 7 is non-negative, where energy is be- 
ing transferred from dark energy into matter. Con- 
cordance models with this modification are shown in 
Figures [7| and m 

In these models, the energy density in dark matter 
peaks before decreasing and eventually becomes neg- 
ative, at which point the models were stopped. The 
density of a component crossing through zero violates 
either the Big Bang condition or the continuity equa- 
tion (|Lindcr 2004) and therefore a simple constant 7, 
while it may give insight into the qualitative effects of 
energy transfer, cannot model an interaction suitably 
over the full range of cosmic time. 

4.2 Dark Energy Decay 



To avoid the problem of negative dark energy density, 
we considered a decaying dark energy model with a 
similiar form as in lTurnen ijlOSiTT) . Here 7 has the form: 



7 = Xpj. 



(29) 



where A is a decay constant. This situation is anal- 
ogous to radioactive decay, with dark energy decaying 
into matter. This model is shown in figures 1^ and llOl 

In these models, the dark energy density peaks be- 
fore decreasing asymptotically to zero. Thus, dark en- 
ergy only dominates the energy density of the universe 
for a finite period of time, and in the future will be 
negligible, just as it was in the past. To illustrate this 
point, we can calculate the deceleration parameter q{t) 
for these models. The deceleration parameter is given 
by: 



lit) 



R\ f R 

R) \R 
l\ fp + 3p 
P 



(30) 
(31) 



by equations 112al and Ill2cll with k = 0. q{t) is de- 
fined so that a decelerating universe has q > 0. An 
Einstein-de-Sitter universe has q = 1/2, a cosmologi- 
cal constant dominated universe has g = — 1, and, in 
general, a universe dominated by an energy component 
with equation of state w has q = (1 + 3«;)/2. 




1(0 =0.00 pA-(0 

■y(t) = 0.25 px{t) 
-(f) =(l.r,|j ,,_A(f) 
7(*)= 0.7.5 px(t) 
7(i) = 1-00 px(i)- 



Time(l/H„) 

Figure 11: The deceleration parameter for a de- 
caying dark energy model: The light blue model 
is the concordance model, where the universe has 
started accelerating in the recent past and will al- 
ways accelerate in the future. By contrast, the 
models with decaying dark energy only experience 
a finite period of accelerated expansion, beginning 
and ending like an Einstein-de-Sitter universe. 



A plot of q{t) is shown in figure [TTI In the early 
universe, matter dominates and the universe is decel- 
erating. The energy density of matter quickly drops 
as the universe expands, and the dominance of dark 
energy causes the universe to accelerate. However, the 
decay of dark energy into matter eventually sees mat- 
ter dominate and the universe decelerates again. Such 
a model would begin and end like an Einstein-de-Sitter 
universe, but would experience an epoch of accelerated 
expansion. 

4.3 The Thermodynamics of Dark 
Energy 

This paper has approached dark energy phenomonel- 
ogically, parametrising its effects on the expansion his- 
tory of the universe without worrying about their phys- 
ical interpretation. However, the ultimate goal of the 
study of dark energy is physical understanding, where 
dark energy models are motivated by a deeper under- 
standing of particle physics. Of particular interest is 
the physical meaning of the interaction parameter 7. 

To this end, consider the adiabatic equation Ijl2bfl . 
It can can be derived from the first law of thermody- 
namics for a closed, constant entropy system: 



dE = d{pV) = -pdV 



(32) 



In our interacting models, we assumed that energy 
was being exchanged between components. Thus we 
consider the first law of thermodynamics as it applies 
to two open systems of Ni particles (or quanta), which 
are exchanging particles. We introduce the chemical 
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potential fi: 



dEi ^-pidV + /lidNi 



(33) 



Comparing this equation to equations 1)27^ and 1128^ , 
we find that our phenomenological parameter 7 can be 
written as: 



i?3 dt 



fj.x dNx 

i?3 



(34) 



As would be expected, 7 is proportional to the rate 
of exchange of particles per unit volume, and can be 
considered the rate of exchange of energy density. 

We now consider the relationship between 7 and 
the dark energy equation of state. By taking gamma 
over to the left hand side, equations 12711 and 12811 be- 



px + 3Hpx{l + wx + 



SHpx' 



= 
= 



(35) 
(36) 



119731 ). Other approache s include alt ering Einstein's 
equations to include f**" , iYilmaz| l982D. sp arking refu- 
tation and counter-refutation llMisnei*199f^:' Alle y et alJ 
.1995,1. and using path groups to formulate Gaussian 
fiux integrals in curved spacetime lMenskv^^2004^ . 

In light of this we must be careful in applying the 
first law of thermodynamics in a cosmological context. 
We previously assumed that energy leaving the dark 
energy component will appear in the matter compo- 
nent. This is equivalent to saying that the universe as 
a whole is a closed system: 



This suggests that we can simply consider interac- 
tion as a particular example of evolving equations of 
state for dark energy and matter. That is, we can de- 
fine the effective dark energy equation of state to be: 
Wx,eB = wx + y/3Hpx and the effective equation of 
state of matter to be: it)m,cff ~ /3H pm. 

But is it a good idea to consider interacting compo- 
nent models to be a special case of evolving equations 
of state? Phenomenologically, the answer is yes. In 
studying the expansion history of the universe, this 
strategy will allow us to examine both kinds of models 
in a common parameter space, knowing that it has no 
effect on the solution of the Friedmann equations for 
a{t). 

However, from the standpoint of physical under- 
standing the answer is no. There is a marked differ- 
ence between physical processes that create matter and 
those that change its equation of state. Matter with a 
non-zero equation of state has pressure comparable to 
its energy density. Pressure is related to momentum, 
so matter with a non-zero equation of state is matter 
moving at relativistic speeds^ . Thus reducing interact- 
ing components to a special case of evolving equations 
of state blurs the distinction between creating pres- 
sureless matter and accelerating matter to relativistic 
speeds. While this has no effect on expansion history 
(a(t)), it will undoubtedly affect structure history, i.e. 
the growth of inhomogeneities. 

One final point can be made from thermodynamic 
considerations of dark energy. Conservation of energy 
in General Relativity is not well understood. Some 
are happy to sa y that energy is simply not conserved 
llHarrisonlll99l) : Einstein introduced the gravitational 
stress energy tensor (t^^) to represent the energy of 
the gravitational field to rest ore energy conservation 
(lEinstein fc GrossmannllTgigf) . but this leads to prob- 
lems because it is a pseudo- tensor: non-localizable, 
non-unique and non-covariant iMisner. Thorne. fc Wheelerl 



^Recall = (pc)^ -i-{mc^)^ . The (mc)-^ term is so large 
for matter that only matter with d — > c has non-negligible 
momentum. 



^diV = M„diV„ + pxdNx = 



(37) 



This seems very reasonable, but the same does not 
apply to the work term {pdV ^ 0). The work done by 
pressure in an expanding universe does not reappear 
in any other form, which is generally interpreted as 
non-conservation of energy in an expanding universe 
(Harrison 1995). Thus, we should keep in mind that 
it may be possible to construct plausible dark energy 
models where H371 does not hold. 



5 Cosmological Implications 

We have discussed cosmological models that explore 
alternatives to a cosmological constant in explaining 
the accelerated expansion of the universe. The testing 
of dark energy models by observational data is a crucial 
challenge for cosmology in the coming years. We now 
outline some consequences of the range of presently 
allowed models. 



5.1 Age of the Universe 

Dynamical dark energy models will give a different ex- 
pansion history R{t) than the concordance model for 
the same present day values of the density parameters. 
While constraints from analysis of the CMB and other 
measures, such as g alaxy distributions, need to be care- 
fuUy considered fsee lOlivares. Atrio-Barandela. fc PavonI 
I2OO5I ■ dynamical dark energy may indicate that the 
age of the universe is different than expected from 
the concordance model. This highlights what should 
be obvious: the age of the universe is model depen- 
de nt. Car e must be taken in this area, for instance 
in iBespr'ov ansv (2005), the parameters of a dark en- 
ergy decay model are constrained by assuming the 
age of the universe which is found using the concor- 
dance model, a different model from that being con- 
strained. Self-consistency in treatment of dark energy 
cosmologies is essential. This demonstrates the way 
in which the figure of 13.7 ± 0.2 Gyrs has become 
broadly applied: this figure is often quoted without 
acknowledging the assumption of a cosmological con- 
stant. To illustrate the increased uncertainty in the 
age of the universe when this assumption is relaxed, 
we find, using the zeroth order constraint on the dark 
en ergy EOS, wx of ~l-24 < wx < —0.74 as found 
in lWang fc Mukherle^ J2004D . that the derived age of 
the universe lies between 12.8 and 13.9 Gyrs, using 
h — 0.7, Wm = 0.3, wx ~ 0.7 as throughout this paper. 
This is not a rigorous analysis, as degeneracies between 
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wx and Q.m are likely to narrow this constraint, but 
does show the order of the increased uncertainty. Fu- 
ture work will determine a more precise value. 

5.2 Coincidence Problem 

It has been proposed that allowing a transfer of en- 
ergy from dark energy to dark matter may hel p re- 
solve the coincidence probl em (see Cai & Wang 20051: 
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A Numerical Model 

The numerical model used to generate the results of 
this paper consisted of a system of first order, coupled 
ODE's which were then integrated in a Runge-Kutta 
scheme. The solution of the Friedmann equations can 
be tackled in two ways: 

We can solve for R in terms of t by solving the 
equations: 



2 SttG, , k 

H = -^(Pm + px) - 

Pm = -3Hpm + 7 
px = -3ii'px(l + wx{z)) - 7 



(38) 

(39) 
(40) 

The alternative is to solve for the energy density 
of the components in terms of z, using the equations: 
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dpm _ 3pm 

dz ^ 1 + 2 ^ (l + z)H{z) 



dpx _ 3pxil +wx{z)) 



dz 



H{z) 



l + z 
8ttG 



+ 



{p,n + px) 



{l + z)H{z) 
k{l + zf 



(41) 
(42) 
(43) 



Then, t as function of z could be obtained by integrat- 
ing: 

dz 



The present is z = 0, and modelling into the past 
is done by allowing z — > oo, and modelling into the 
future is done by allowing z —1. R{t) can then be 
recovered from z{t). 

This second approach makes it easier to calculate 
the angular diameter and related distances. This could 
be done using the equations in section |5| However, 
the equations for angular diameter distance have been 
developed in the case where the universe is "clumpy" 
i.e. where it is only homogeneous and isotropic on 
average. The resultant equation is the Dyer-Roeder 
equation, as generalised in iLindcr ( 1988.'! : 



s + q{z) 

l + z 



rd + 



(^(TTI) + 3^..(2))a,(2)n,(2: 



) Ud = (45) 



where q{z) is the deceleration parameter, and the ini- 
tial conditions are: 



rd{zo, zo) = 
dr(zo, z) I _ Ho 



dz 



H{zo){l + zo) 



(46) 
(47) 



where we have deliberately not set zo = so that an- 
gular diameter distance can be calculated for observers 
at any redshift, as is common in situations involving 
gravitational lensing. The dumpiness parameter (ai), 
defined to be the ratio of the amount of component i 
smoothly distributed to the total. It is set to 1 in this 
paper, which gives the same solutions for angular di- 
ameter distance as section |21 Our codes can relax this 
condition, and future papers may explore the effect of 
inhomogeneities in dark energy models. 



